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ABSTRACT 

We analytically determine the correlation functions of the stochastic 
q\ , response of a generic mapping system driven by colored noise. We also ad- 

a: 

dress the issue of noise cascading in coupled-element systems, particularly 

^ . 

^ , in a uni-directionally coupled system with a dc-filtered coupling. A noise- 

feedback method and an opposite-pairing method are introduced to control 

O 

the stochastic response. Excellent agreements are obtained between analyti- 
cal results and numerical simulations. 
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1 Introduction 



The stochastic response of dynamical systems to noise has attracted much 
attention recently and is a recurrent theme in the study of nonlinear sys- 
tems. Some classic works were accomplished on the persistence of dynamical 
system under random perturbation Jl|. Recently, there are many novel and 
intriguing phenomena induced by noise, for instance, noise-induced order 
, synchronization by external noise || , stochastic resonance , noise-driven 
coherence resonance || , and noise-assisted stabilization |J. Some experi- 
mental verifications of these phenomena have been found in extremely com- 
plicated real systems, such as the brain [7J] and plasma ||, and the physical 
picture can be captured by dynamical models with simple additive noise. On 
the other hand, noise-induced multistability has been addressed in explor- 
ing the influence of multiplicative and additive noise M, and the reentrance 
phenomenon was introduced in studying the dynamical system with colored 



noise 



101. The colored noise is the noise source with finite correlation time. 



In opposite to colored noise, white noise have not any correlation between 
time lapse. Whenever the noise is dynamical involved, the correlation func- 
tions of the stochastic response of dynamical system to colored noise is usu- 
ally very difficult to obtain exactly. Several theories, for examples, the small 
correlation time theory [|11|], the functional-calculus theory P2| , the decou- 
pling theory |T5| |, the unified colored- noise theory JTJJ and the interpolation 
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procedure JT5|, were proposed to treat the colored noise problem in stochas- 
tic processes. However, analytical results can only be obtained under the 
limit of small or infinite correlation time. Therefore, how to determine the 
stochastic response of the colored noise exactly, is still an important issue. 

In another front, it is known that collective behaviour which arises from 
the coupling among the elements could be dramatically different from the 
behaviour of individual element, and is also a current interest in nonlinear 
dynamics. In the deterministic systems, some novel features, such as self- 
induced spatial disorder and chaotic itinerancy, have been numerically found 
in distributed nonlinear element systems |16|]. Universal critical behaviours 
have also been found in coupled map lattices |I7j and the first order phase 
transition was studied in coupled oscillator systems JjlJ]. Nevertheless, to- 
ward a more realistic model, one should include the noisy elements. Indeed, 
the stochastic resonance has been explored by numerical simulation and ex- 
periment in coupled-element systems with additive noise JTDJ. However, it 
is of fundamental interest to provide more illustrative example to clarify the 
exact role played by noise in the coupled system. 

In this paper, we analytically determine the correlation functions of the 
stochastic response of the steady state when the colored noise is dynamically 
involved in a generic mapping system. We also investigate the issue of noise 
cascading in coupled dynamical systems. As an analytical illustration, we 
adopt a special model in which all the elements are uni-directionlly coupled 
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to each other, specifically, we consider a dc-filtered coupling. The coupling is 
fundamental different to the global or nearest-neighbour coupling in spatially 
extended system[16,17], which is receiving much attention nowadays. But, 
this simple model thus serves as an analytically solvable model to illustrate 
the nature of noise during propagation in spatial extent under the influence 
of dynamical structures. Furthermore, we found two elegant ways to control 
the stochastic response. 

2 Stochastic response of colored noise 

To begin with, we consider a generic one-dimension nonlinear map in 
which the noise is involved dynamically, 

x n+l = fr( x n) + Cn- (1) 

where ( n can be colored or white noise with zero mean fluctuation. For 
simplicity, we only treat the stochastic response of the fixed point, one can 
extend our results to higher periods by following the similar procedure in 



Reference |f20| . The dispersive output can also be expressed as 

x n = x 4~ £ ra , (2) 

where x is the noise-free fixed point, x = f r (x), and £ n is the stochastic 
response. From now, the amplitude of the noise and stochastic response 
are assumed to be small in comparison with the fixed point, such that the 



4 



perturbation scheme is applicable. By substituting Eq.(2) into Eq.(l), a 
recurrence relation is obtained: 

= M^n + Cn, (3) 

where M = ^r%- With \M\ < 1, we will , after the transient, have a small 
response fluctuation around the stable noise-free fixed point. By iterating 
Eq.(3), we have 

k 

& = M fc £o + E W -1 c*-i (4) 

After passing the transient state, i.e. k » 1, the first term of Eq.(4) can be 
neglected. As a result, the mean fluctuation of stochastic response is zero, 
if < (k >= 0. Here, < A k >= limAr^oo-^X^Q 1 ^Wi denotes the long time 
average of A k+i . The index k just reminds us that we are taking average for 
those states after the transient k » 1, and it is irrelevant to the average. 
The zero mean of stochastic response imply that the average of dispersive 
output x n is same as the fixed point x. The general relations between the 
correlations of stochastic response and those of input noise can be obtained 
in compact forms: 

< tek >= {< CfeCfc >+2E M '< (k(k + i > j (5) 

< z k £ k+t >= X _ M2 1< CfcCk >+2E F< CkC k+l >| 

l oo 

+EE^ 3+i <aw ! >, t>i. (6) 

j=l i=0 
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If the correlations of the input noise are given, one can obtain the correlation 
of stochastic response after the summations are accomplished. For example, 
if the input noise £ fc is a white noise and is denoted by r/ k which have zero 
mean and derivation D =< rf >, then the mean fluctuation of stochastic 
response is zero and the correlation function of the stochastic response is 

M e D 



£>0. 



(7) 



1 — M 2 ' 

Eq.(7) shows that the stochastic responses are colored when the white noise 
is dynamically convoluted by a generic map unless the system is critical i.e. 
M — 0. The correlation time r c of the colored response can be identified as 
r c = -(|ln|M||)- 1 . 

Next, we use the result of Eq.(7) as the input colored noise to drive an- 
other system, x n+ \ = f r (x n )+^ n . With Eqs.(5) and (6), the mean fluctuation 
and the correlation of the stochastic responses, £ n = x n — x, of the new system 
could be deduced. They are < >= 0, and 



< > = 



< ikik+i > — 



D 



[1 + MM) 



l-M 2 ) [(l-M 2 ) (l-MM)J ' 

d r i 



(8) 



x 



l-M 2 ) 1(1 — MM) 

' t M l — M l -,,1 + MM, 
M(-^ — ) + M\- ) 



M-M 



1 - M 2 



£>1, (9) 



where \M\ = < 1 and x = f r (x). The response and the correlation 

time are related to the dynamical structures M and M explicitly. Eqs.(8) 
and (9) also suggest that the correlation of the response could be greatly 
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enhanced by the variation of dynamical structure. Referring to Fig.l, when 
the first stage is at the critical point (M = 0), the stochastic response of 
the second stage exhibits a symmetry at the critical point (M = 0). The 
symmetry is broken as the dynamics of the first stage moves away from its 
critical point. The minimum of stochastic response will shift to MM < 
region. This fact will provide us a way to control the stochastic response and 
we will return to this point later. 

Moreover, The extention of this result to higher periods case straightfor- 
ward but tedious. For example, to explore the stochastic response of each 
state of period-2, we shall divide the out data into two sets, called even and 
odd parts, and estimate the correlation for each set individually. Then, the 
correlations of stochastic response on period-2 states can be simply carried 
out by converting M and ( n to M 2 = §£ \ Sa and b n = |f | 5 < n + Cn+i in 
Eqs.(5) and (6) respectively. Where x 1 ( x 2 ) is the mean value of odd (even) 
part of dispersive output, i.e. the period-2 orbit for noise-free case, and x' is 
x 2 ( xi) when the index n of b n is even (odd). 

3 Noise cascading 

As mentioned above, the noise cascading problem in coupled-element 
systems, especially for spatially extend coupled system, is remained as a 
challenge issue to be addressed. To gain more illustrative features, we take a 
dynamical system, called the first stage, which is driven by a white noise. The 
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stochastic response is extracted by a de-filter and used as the noise to drive 
another system, called the second stage. This procedure could be repeated 
successively. Physically all the dynamics are described by the following 
mapping, 

r (1) fQ-UrW} +r? 

x n+l — Jr \ x n J ' 'In 

41 = fP(^)+e (io) 

M f(m)( (m)\ , Am-1) 

x n+l ~~ Jr l x n ) ' Sn ) 

where = x$ — , and can be extracted by a dc- filter. The essential 
feature of noise cascading can be addressed by treating identical systems (all 
are the same). To obtain an explicit form of analytical result we will 
assume \M\ « 1. Using Eqs. (5) and (6), and keeping the lowest order 
of M, it can be proven that the correlation functions for the m-th (m > 2) 
stage are 

<d m) d m) > - (l + m 2 M*)D (11) 
<^ m) £K> « C%r*M'D, t>l, (12) 

where C^t!^ -1 is the binomial coefficient. The proof is shown in appendix. 
These results indicate that the standard deviation of the noise is amplified 
when more stages are included. As shown in Fig.2, the analytical result 
is also supported by numerical simulation . We note that the correlation 
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functions monotonously decrease with i for M > and oscillating decay for 
M < 0. Both damping factors have r c = -(| In \M\ for £ » I. 

4 Control of noise response 

The noise amplification found above eventually leads to the breakdown of 
the coupled-element system as long as there are too many elements . This 
amplification phenomenon raises an important issue, namely, whether or not 
a controllable stochastic response can he obtained in noise cascading. If the 
answer is negative then it is hopeless to study noise cascading in dynamical 
systems. Fortunately, we found two way to surmount the above difficulty 
The first solution to this problem can be achieved by means of feedback. 
Let us modify the generic one-dimension map, Eq.(l), by adding feedback 
response a£ n , 



where a is the feedback strength and £ n = y n —y. By using the same approach 
as before, we find the modified general relation between the correlations of 
stochastic response and those of input noise, 



Un+l = friVn) + (n + <, 



(13) 



< CkCk > 



1 — (M + a) 



1 — (M + a) 
(M + a) e 



1 




£ oo 



+ ^^(M + ap ! <( fc a +J+! >, i>\. (15) 



j=\ i=0 
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This result indicates that the output response at best could be reduced to the 
input noise if we adjust a to be o = —M. Fig. 3 is the simulation result and 
shows that the feedback scheme can reduce the stochastic response efficiently. 
The standard derivation of output response is minimized at a = —M, and 
it is indeed smaller than the case without feedback (a = 0). On the other 
hand, as implied by the results of Eqs.(8) and (9), some suitable variations 
of dynamical structures may be employed to reduce the stochastic response, 
and thus postpone the breakdown. For a simple illustration, let us consider a 
two-stage system with opposite pairs. It means that we choose the stability 
quantity of second stage to be M = —M, which is opposite to the previous 
one. By using Eq.(8) and (9), the correlation of the output noise is 

M 2e D 

and < £,k£,k+2e+i >= 0. Since \M\ < 1, it is clear that the output noise is 
suppressed in comparison with Eq.(7) and the correlation time is still r c = 
— (|ln \M\ in this case. That offers us an opportunity to reduce the colored 
noise by choosing the dynamical structure with opposite pairs. Extending 
this opposite-pairing method to the cascading model with p-pair, M ( 2p_1 ) = 
M and M<&) = —M, p > 1, one finds that the correlation functions for the 
p-th pair are 

<e?" 1) Cf" 1) > « (1 + M 2 +/M 4 )D (17) 
<d 2p) d 2p) > « (1+P 2 M*)D. (18) 
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These results can be proven by following the same procedure shown in the 
appendix. As shown in Fig.4, the output noise is indeed greatly reduced and 
the breakdown is postponed. 

5 Concluding remarks 

In summary, we analytically derive the relation between the correlation of 
stochastic response of a steady state in a mapping system and that of input 
noise. We also obtain the exact form of the correlations of the noise cascading 
in an uni-directionally dc-filtered coupled system. Our results suggest that 
noise amplification might lead to catastrophic disaster in coupled-element 
systems. However, in certain cases, such as the one discussed above, the 
amplification can be controlled either by a feedback scheme or by an opposite- 
pairing arrangement of dynamical structures. Let us address the implication 
on this point. Assume that the noise-cascading system will breakdown when 
the derivation of stochastic response is large than a threshold D t h, thus the 
noise amplification set the limit on the cascading length. For examples, the 
maximum lengths of the identical-coupled and opposite-pairing cascading 
models are m max = [\j {Dth D D) pgy], and 2p max = 2[^ (Dth ~ D) jjp] respectively, 
where [ ] is the Gaussian notation. Since \M\ < 1, we will have a shorter 
chain for identical elements and a longer chain for the opposite-pairing one. 
On the other hand, the noise-feedback coupling will give an un-limited chain. 
Thus, it is worthwhile noting that our models, even they are not so realistic 
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now, may provide a scheme in which the noise could influence the size of the 



pattern formation pi 



Finally, there are some remarks should be mentioned. All the analytical 
results agree with the numerical simulations very well as long as the ratio 
between the amplitude of input noise and the value of fixed point is less than 
5%, and the dynamical systems are not too close to the bifurcation point. 
The same conclusion can also be drawn for the Ikeda map. It is important to 
note that the dc-filtered coupling is essential through out this investigation. 
Without the filter, the dynamical structures of each stage will be greatly 
altered by the dc-signal and the coupled system will quickly breakdown. 
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Appendix 



Here, we will use mathematical induction to prove Eqs.(ll) and (12). 
After appling the approximation \M\ « 1 and keeping those leading terms, 
Eqs. (5)- (6) end up with two general relations between the (m + l)-th stage 
correlation functions and m-th stage correlation functions, i.e., 

< d m) ei m) > « (1 + M 2 ) < ef- 1 ^ > +2M < it n ' 1] & >, (19) 

+E^ ! <?r , d:- l) >, €>i. (20) 
i=i 

We note that the leading order of < ^'^Zk+i^ > * s ; . 

Let us first consider the m = 2 case. We read out < ^jp^ >~ (1 + 
M 2 )M l D from Eq.(7), and put it into Eqs.(19) and (20). It is easy to check 
that 

<tftf> « (1+4M 2 )D (21) 
< ^k+e > ~ ^ + 1M ^ D = Cf +1 M'D., £ > 1. (22) 

Next, for the cases of m > 2, we will show that if Eqs. (11) and (12) are 
true for m-th stage then they are also true for the m + 1 case. Plugging Eqs. 
(11)-(12) into Eqs. (19)-(20) and expanding them to those leading order 
terms, we obtain 

< ef +1) ei m+1) > ~ (1 + M 2 )(l + m 2 M 2 )D + 2C™_ 1 M 2 D 

13 



« (1 + (m + l) 2 M 2 )D (23) 
<ef +1) & +1) > ~ M i {l + m 2 M 2 )D + Y j C i +X 1 M"D 



i=0 



i=i 

E <#r ^ = c^+t^A ^ > i, (24) 



where the summation identity of binomial coefficients is used. By mathe- 
matical induction, Eqs. (11) and (12) are true for any m-th stage. Q.E.D. 
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Figure Captions 



Figure 1. 

Up: The block diagram for the stochastic response £ n of the dynamical sys- 
tem M where the correlation function of the input noise is Eq.(7). That is 
equivalent to a uni-directionally coupled system in which two elements are 
coupled by a dc-filter and with white noise input. 

/ < £2 > 

Bottom: y via M for several fixed M. The solid lines denote the an- 
alytical estimations. The circle (bullet, box) denote the simulation results 
for the dynamical structures M = 0.5 (0.0, —0.5) of first stage. Here, the 
dynamical systems are logistic maps f r (x) = rx(l — x) with M = 2 — r, and 
fr(x) = fx(l — x) with M — 2 — f. The input uniform white noise r/ n have 
amplitude 0.01. 
Figure 2. 

Up: The block diagram for an identical coupled-element system. 

/<£ (m)2 > 

Bottom: y — ^ — via stage number m, for two identical coupled-element 
systems which have different elements f r (x) = 1.9 * x(l — x) and f r (x) = 
1.99 * x(l — x), respectively. The solid lines denote the analytical results 
Eq.(ll). The bullet ( box) denote the simulation results for the input noises 
with amplitudes 0.01 (0.001). 
Figure 3. 

Up: The block diagram for a feedback scheme to reduce the stochastic re- 
sponse. 
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Bottom: y via the feedback strength a. The input uniform white 
noise has amplitude 0.01. Two logistic maps f r (x) = 1.8 * x(l — x) with 
M = 0.2 and f r (x) = 2.3 * x(l — x) with M = —0.3 are demonstrated. Both 
simulations and analytical results show that the responses are minimized at 
a = —0.2 and a = 0.3 respectively. 
Figure 4. 

Up: The block diagram for a opposite-pairing scheme to reduce the stochas- 
tic response. 

/<£ (m)2 > 

Bottom: y — ^ — via stage number m, for the opposite-pairing coupled 
system which have a pair of logistic maps f r (x) = 1.9 * x(l — x) and f r (x) = 
2.1*x(l — x). The solid lines denote the analytical results Eqs.(17) and (18). 
The bullet denote the simulation result for the input noise with amplitudes 
0.01 . 
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